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ABSTRACT 

All monochromatic opacity data from the Opacity Project (OP), together with all 
codes required for the calculation of mean opacities and radiative accelerations for 
any required chemical mixture, temperature and mass-density, are being put on a 700 
Mb CD which will be made generally available. The present paper gives a concise 
summary of the contents of the CD. More complete documentation will be provided 
on the CD itself. 
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1 INTRODUCTION 

Methods used by the Opacity Project (OP) for the calcu- 
lation of mean opacities (Planck and Rosseland) were de- 
scribed in Paper I (Seaton et al., 1994); further work on the 
inclusion of contributions from inner-shell transitions were 
discussed in Paper II (Badnell and Seaton, 2003) for the 
six-element mix of Iglesias and Rogers (1995); comparisons 
with OPAL (see Rogers & Iglesias, 1992) for that mix in 
Paper III (Seaton and Badnell, 2004); and up-dated results 
for all elements, with inclusion of inner-shell contributions, 
in Paper IV (Badnell et al, 2004). The OP approach to the 
calculation of radiative accelerations is discussed by Seaton 
(1997). 

The basic monochromatic opacity data for the 17 el- 
ements H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, 
Cr, Mn, Fe and Ni occupy 665 Mb as gzip formatted files, 
or 868 Mb as binary files. Some improvements have been 
made in the frequency-mesh which is used (see section 2.2.1 
below) which lead to improved accuracies in the numerical 
integrations required to obtain mean opacities and radiative 
accelerations. It is found that the integrations for Rosseland- 
means never give numerical errors larger than 0.1%. Im- 
provements have been made in the codes used to obtain 
mean opacities and accelerations for any required chemical 
mixture, interpolated to any required values of temperature 
and mass-density. The interpolated Rosseland mean opaci- 
ties are shown to be numerically correct to better than 1% 
for all practical purposes. 

It is our intention to make all OP data generally avail- 
able, both on a 700 Mb CD which will be described in the 
present paper and on a web OPserver which will be described 
by Mendoza et al. (2004). 

The CD contains all data required for the calculation 
of OP mean opacities and radiative accelerations, together 



with full documentation and all required codes, and is avail- 
able on request 1 . It contains a tar file of compressed (gzip) 
files. The file OPCD.ps.gz provides full documentation, in- 
cluding information about processing of all other data on 
the tar file. 

The user is urged to read the file 0PCD .ps before decom- 
pressing any other files. 

The present paper gives a concise summary of the con- 
tents of the CD. Full documentation is given in the file 
0PCD.ps. The CD includes input and output files for some 
test runs. 



2 BASIC DEFINITIONS 

2.1 Monochromatic opacities 

The monochromatic opacity cross-section for element k is 
ak(u) where u = hv/iksT) and fee is Boltzmann's constant. 
The correction-factor for stimulated emission, [1 — exp(— u)], 
is not included 2 in the definition of a^. For a mixture of 
elements let element k have abundance by number-fraction 
of /fe, fk = 1- The cross-section for the mixture is 

o-(u) = 22 fk°~k(u). (1) 

fc 

2.2 The Rosseland mean 

The Rosseland-mean cross-section is 

1 Requests to: Claude Zeippen (Claude.Zeippen@obspm.fr); or 
Anil Pradhan( pradhan.l@osu.edu) 

2 We here consider the quantity actually tabulated in the OP 
archives. A number of previous papers have defined <r fe with in- 
clusion of the stimulated emission factor. 
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where 



F(u) = [15/(47r 4 )]it 4 exp(-u)/[l - cxp(-u)] 2 . 



(2) 



(3) 



The Rosseland-mean opacity, per unit mass, is kr = an/n 
where /i is mean atomic weight. 

Numerical values for cross-sections a will be given in 
atomic units (a 2 ) and values for opacities kr in cgs units 
(cm 2 g- 1 ). 



2.2.1 The frequency variable 
It is convenient to use the frequency- variable 3 
F(u) 



v(u) 



giving 



o [1 - exp(-w)] 
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(5) 



where Voo = w(w — > oo) 



2.3 Radiative accelerations 

The radiative acceleration for element k in a stellar interior 
is 



5rad = (l/c)(n/nk) Fun-lk 



(6) 



where c is the speed of light, jik the atomic mass of element 
k and 



F = nB(T eS )(R i ,/rf 

where T e a is the effective temperature, 

B{T) = 2(nk B T) 4 /(l5c 2 h 3 ), 



(7) 



(8) 



is the radius of the star and r the distance from the 
centre of the star. In (6), 
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(9) 



where cr™ a is a cross-section for momentum-transfer to the 
atom, 



cr™ ta = <Tk{u)[\ - exp(-tt)] - afc(u). 



(10) 



In (10), dk{u) is a correction factor to take out contribu- 
tions from scattering processes and momentum transfer to 
electrons in photo-ionization processes. 

Values of g r ad are given in cgs units (cm s~ 2 ). The quan- 
tity 7fe is dimensionless. 

One is usually interested in the case in which the abun- 
dance of the selected element k varies with depth in a star 
due to the operation of diffusion processes. The codes allow 
for the possibility of varying the abundance of k by a fac- 
tor x, keeping the relative abundances of all other elements 
unchanged. 

3 Some previous papers have taken a to be defined with inclusion 
of the stimulated emission factor, and v to be defined by v = 
J F(u) du. 



3 ARCHIVED DATA FILES 

All data are tabulated on a mesh of values of temperature T 
and electron-density N e . The indices ite and jne are defined 
by 



ite = 40 x log(T), jne = 4 x log(iV e 



(11) 



All data on the CD are with the "m" mesh, even values of 
ite and jne, Aite=Ajne= 2. 

The cross-sections are tabulated for ntot= 10000 
equally-spaced values of v. The <Jk data are in files mzz.ttt 
and the are in files azz . ttt where zz is a two-digit num- 
ber specifying nuclear charge (zz= 01 for H) and ttt is 
a three-digit number specifying ite. Data for all jne in- 
cluded are contained in each files mzz.ttt or azz. ttt. Files 
mzz.mesh give values of u for the tabulated values of v. 

The OP codes use unformatted (binary) data files, for 
economy in storage and speed of input/output. However, 
the CD contains files which are formatted and compressed 
(gzip) for ease of transfer between machines. The proce- 
dures to be used for converting to binary are described in 
the file OPCD.ps. 

Some summary data, described in the documentation, 
are given in files mzz. index, mzz.ion, mzz.smry. 



4 CODES FOR ROSSELAND-MEANS 
4.1 The code mixv.f 

The code mixv . f reads a list of abundances fk for a mixture, 
obtains the cross-sections a and produces an output file with 
values of mass-density p and Rosseland means kr for all 
(T, N e ) mesh points. 



4.2 The code opfit.f 

The code opf it . f reads an output file from mixv . f and a list 
of values of (log(T), log(p)) for a stellar model and produces 
a file giving values of 



log(K R ), 91og(/tR.)/51og(T), 91og(/tR)/51og(p). 



(12) 



The code uses refined techniques of bi-cubic spline interpola- 
tions. It has an optional facility (rarely required) to smooth 
the data before doing the interpolations, and an optional 
feature to produce output tables in OPAL format. 



4.3 The code mx.f 

The code mx . f produces output similar to that from opf it . f 
in a one-step process, using less refined bicubic interpolation 
techniques. It has a facility that abundances can be changed 
at any point during a run for a particular stellar model. 



4.4 The code mixz.f 

The code mixz . f produces monochromatic opacities for any 
specified mixture of metals. The output files are stored as 
a "fictitious" metal, mzz.ttt with zz a two-digit number 
greater than 28 (data for Nickel are stored in files m28.ttt). 



OP data on CD for mean opacities and radiative accelerations 3 



4.5 The code mxz.f 

The code mxz.f is similar to mx.f but uses the files for a 
metal-mixture from mixz.f. It has a facility that, at each 
depth point, the values of X and Z can be changed. 

4.6 Choice of mixture codes 

Use of mixv.f plus opfit.f can be convenient if the same 
mixture is to be used for many stellar models. The use of 
splines in opfit.f ensures that the derivatives in (12) are 
continuous. It is shown in Paper IV that, for all cases likely 
to be of practical interest, the interpolations using opfit.f 
do not introduce any errors as large as 1%. 

Use of mx . f is particularly convenient if there are 
changes in composition depending on depth in a stellar 
model. In occasional difficult cases use of mx . f can give a 
lack of smoothness in the derivatives. 

The test runs of mixv.f, opfit.f and mx.f are for 
S92 abundances (see Paper I), log(T) = 6.0 to 7.3 and 
log(.R) = —1.75, roughly corresponding to conditions in the 
solar radiative interior. Use of mx.f gives values of log(/-CR.) 
and <91og(ftR,)/<91og(T) which appear to be smooth and in 
agreement with values from opfit.f but there is some ap- 
parent lack of smoothness in <91og(ftR,)/i91og(p) from mx.f, 
as shown on Figure 1. 

The codes mixz . f and mxz . f were used in Paper IV for 
the case a more realistic solar model in which, below the 
base of the convection zone, X and Z varied continuously 
as functions of depth. 



5.4 Choice of codes for accelerations 

Use of accfit.f provides the most accurate interpolations. 
The code ax . f can be useful for cases in which abundances 
vary with depth. 

The test case is for Argon in a stellar model with T c g = 
10 4 K. Figure 2 shows p rad for \ = 0.01, 1.00 and 100. The 
deep minimum in the vicinity of log(T) = 5.45 occurs when 
the Argon is going through a Neon-like ionization stage. It 
is seen that g ra d from ax.f shows some lack of smoothness 
for the rather extreme case of x = 0.01 
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5 CODES FOR RADIATIVE ACCELERATIONS 
5.1 The code accv.f 

The code accv . f has a structure similar to that of mixv . f . It 
reads a list of abundances and specification of the selected 
element k, reads data specifying a range of values of the 
abundance-multiplier x, reads the files mzz . ttt and azz . ttt 
and produces an output file giving, for each (T,N e ) mesh 
point, values of 

kr, dun/dx, 7, <h/dx (13) 
for each value of x- 



5.2 The code accfit.f 

The code accfit.f is similar to opfit.f using bi-cubic 
splines. It reads a file output from accv.f, a list of values 
of x, the effective temperature T e g for a stellar model and 
a list of values of log(T), log(p) and r/R+. It produces a file 
giving value of 

log(KR), log( 7 ), log(p rad ) (14) 
as functions of x- 



5.3 The code ax.f 

The code ax.f has a structure similar to that of mx.f. Its 
use is convenient when the mixture of elements depends on 
depth in a stellar model. 
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Figure 1 <91og(KR,)/i91og( / o) for a model with S92 abun- 
dances and log(-R) = —1.75. Results obtained using the 
codes opfit.f and mx.f. 
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Figure 2 Radiative accelerations for Ar in a model with 
T e ff = 10 4 K. Results from codes accv . f and ax . f with abun- 
dance multipliers x = 0.01, 1.0 and 100. 



